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ABSTRACT 

A single dimensional sinusoidal parameter perturbation adaptive 
control system is studied. The operation and signals at every point of 
the system are analyzed theoretically. Then, the system is simulated on 
analog computer starting with the simplest circuit and adding more com- 
ponents until a complete system is formed. The results of simulation 
study are compared with the theoretical analysis. Both agree in most 
cases, Finally, an analytical design of the system is attempted based 
on the simulation study and using a linearizing technique. It was simu- 
lated on analog computer and the performance was compared with the one 
predicted analytically. It does not follow exactly as predicted but the 
deviation is not so great. This method is useful in obtaining preliminary 
design which will be refined by experimental verifications, 

The writer wishes to express his appreciation for the assistance and 
encouragement given him by Professor К. С. Dorf of the U. S. Naval Post- 


graduate School in this investigation. 
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1. Introduction 

The adaptive control system has received much attention recently. 

The principle of adaptive control 13 to maintain a system always in an 
optimum operating condition whatever the affecting situation may be. There 
are two kinds of adaptive control systems; the signal adaptive system and 
the system adaptive system. The signal adaptive system is a system which 
can analyze incoming signals and automatically adjust the system to give 

an optimum output. The system adaptive system is a system which can mea- 
sure its own characteristics and then change its structure or design so as 
to continuously maintain a desired performance. This paper deals with only 
the latter, The system adaptive system should have the following two es- 
sential features to achieve its goal; 

(1) The measurement of its own characteristics. 

(2) The generation of proper corrective signal for the control- 
led system according to the measurement. In the measurement of the 
characteristics of the system, the optimum operating condition or figure 
of merit is defined. The figure of merit or index of performance (IP.) is 
an indication of the characteristic or quality of a system. The IP may be 
a physically meaningful value such as heat produced per gallon of fuel or 
theoretical quantity such as RMS error or mean squared error. 

The sinusoidal parameter perturbation adaptive control system is a 
method of system adaptive control and in this paper the operation and 
signals at every point of this system, are analyzed theoretically. Then, 
the system is simulated on analog computer starting from the simplest cir- 
cuit. Adding more components one by one, the effect of each component is 
observed until the complete system is formed. The results of these simula- 


tion study are compared with the theoretical analysis. Both agree in most 
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cases. Finally, a design of the system is attempted analytically based on 
the simulation study and using a linearizing technique. It was also simu- 
lated on the analog and the digital computer and the actual performance of 
the system is compared with the predicted performance from the analytícal 
design. The results show that the actual performance does not exactly agree 
with the performance predicted analytically but it follows reasonably close 
pattern predicted analytically. This is expected because the analytical 
method involves approximation, This method is useful to obtain a prelimin- 
ary design which will be refined by experimental verification and gives a 


direction of improving the first trial design. 





2, The Sinusoidal Parameter Perturbation Adaptive Control System. 

A single dimensional sinusoidal perturbation system is as Fig. 1. 
The model is an optimum operating model. When the plant is operating 
exactly as the model there will be no error output. But when a parameter 
in the plant is changed to some other from the optimum parameter value due 
to environment change, there will be an error output. By squaring the er- 
ror we get a non-negative function of error. As an index of performance 


we may adopt average of error squared as; 


7 
ТР = 2 zd (2-1) 


This quantity is a function of the characteristics of input to the plant 
and also of the characteristic of plant parameter drift caused by environ- 
ment change. Fig. 2 shows the contours of IP(X,Y) for two dimensional 
case where IP is a function of two parameters X and Y. The minimum IP 
occurs at X=X, and Y= 0’ ТЕ we fix the value of Y at Yo and vary only X, 
the IP variation with X becomes as Fig. 3, the single dimensional case. 
Generally the IP with one parameter variable tends to be a parabolic func- 
tion of the parameter and can be approximated near the optimum parameter 


value as; 


IP(X) = a (X-X) +a (252) 


2 


where Xo is optimum parameter value and а, апа و8‎ 
are arbitrary constants which change the shape 
of IP curve such that it represents closely the 
actual IP variation of any plant. 

In some plants a parameter has a slight or no change with respect to the 


environment change and such parameter is fixed at the optimum value and the 
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C.SINGO t 


IP 
PARAMETER PARAMETER MEASURING 
PERTURBATION CORRECTION DEVICE 
INPUT SIGNAL 
MEASURED IP 
DETECTION REFERENCE 
(a) General set-up 
MODEL 
INPUT BAND PASS 
SQUARER FILTER 
+ 
OUTPUT 
C_SINW,t 4 _ SMOOTHING | 
2 2 FILTER DETECTOR 
PHASE LAG 
FILTER 


(b) Model Reference System 


Fig. 1. Single Dimensional Sinusoidal Parameter Perturbation Adaptive 
Control System. 
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Fig. 3. IP v.s. parameter X 





Contours of IP(X,Y) for two dimensional case. 


Fig. 2. 
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number of parameters to be controlled is reduced. 

Referring to Fig. 3, suppose we perturb the parameter X sinusoidally 
by a perturbation signal, С = с, біп wit. If the actual plant para- 
meter X is at the optimum parameter value, Xo? there will be no IP varia- 
tion output at frequency W,. If the X value is lower than Хо, there will 
be an IP variation output at frequency e and the output signal will be 
180° out of phase with respect to the input perturbation signal as seen in 


the Fig., on the contrary, if the X is higher than the X., the output is in 


0? 
phase with the input perturbation signal. It is also noted that the ampli- 
tude of the output is proprotional to |x-x.|, therefore the output is consid- 
ered to be a carrier, modulated by the variation of | X-X,| through the IP 
characteristic. The frequency of carrier is e, and it carries the informa- 
tion about the polarity of the X-X. by its phase difference with the input 
perturbation signal. When this output is demodulated by multipling with 

the original perturbation signal, it will produce a D.C. signal with proper 
sense to correct the plant parameter deviation. 

The bandpass filter eliminates unwanted signal components contained in 
the modulated signal and also affects the response of the adaptive loop to 
a step change-of signal. 

The detector is a multiplier which multiplies the modulated signal with 
original perturbation signal and produces positive or negative correction 
signal. 

The filter after the detector is a smoothing filter which filters out 
high frequency components in the correction signal and acts as compensator 
to stabilize the adaptive loop. 

The integrator is a D.C. motor in practical operation and it integrates 


correction signal and drives the plant parameter control to the optimum value. 
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If we want to make more than one parameter adaptive, separate loops 
for each different parameter must be established and use different perturba- 
tion frequency for each parameter. The same error measure can be used for 
each loop if it is suitable, but for each parameter, other elements must be 
duplicated, | 

Sometimes, a plant has several local minimum in its IP characteristic. 
In that case, a special measure is taken to find the true minimum, But 
generally, all plants have a single minimum around the optimum parameter 
value. The sinusoidal perturbation method can also be applicable to an 


adaptive system which searches for a maximum IP value. 








3. The Analysis of Signals in the Adaptive System 
3.1 Block diagram of an adaptive system 

To set up a block diagram of a single dimensional parameter perturba- 
tion adaptive control system, the circuit models of each component are 
established. Referring to Fig. 1, the parameter perturbation signal per- 
turbs a plant parameter and the parameter itself also is changing its value 
due to environment change and parameter correction signal coming through 
adaptive loop. This parameter perturbation signal and parameter drift 
cause the error output while the plant and model are operating with input, 
К. It is necessary to find out an equivalent transfer function block be- 
tween the parameter perturbation input and the error output. The para- 
meter drift is also an input for the error output and the input is same as 
perturbation input. We may designate this equivalent transfer function as 


error measurement transfer function and express as 


m 
6268) = ТИ в 8 


K_ and СА are decided by the signal transfer characteristic between 
perturbation input and error output. To find out the transfer character- 
istic, phase lag and attenuation introduced to this error output with res- 
pect to the original perturbation and parameter drift are observed. Al- 
though the perturbation and drift do not come into the plant through 
ordinary input point but through parameter itself, these signals encounter 
the same phase shift as a signal of the same frequency coming into the 

plant through the input of the plant, encounteres. This will be confirmed 
by actual simulation study in Section 4.1. The attenuation factor of this 
error signal is not the same as the attenuation factor of the signal of same 


frequency coming through plant input. The attenuation factor to the error 



















-— N 
> — -- “Өне жок өз өш) | 
ж жы- “-- ы. s-r. 
ae A % ез) 1 سوک‎ ж ccm 
u ..-. U. — 

— — t w < $ @ 
——À — c 6-0 ри 
— — —— — 
- oma cac ps - 4 шы ےھ‎ u U 
—- — s ٭ لی‎ ~ 
= m a ШШ 


(- -%-.,- с 


e ea anal رو موس‎ eee دص گے نے‎ Й 











ar & "c чш». جج مه‎ 
— J 1 мәне وس ستتھط۱ن۔‎ фә - ж 
i = of c— = — w d me a 
Г ТЕЕ e. sr hos ¿e ee 
— — — СТЕП СЕСИИ ЕЕЕ 
w ж A-—: 9 9 o нкф wu 
— <= i) ome l — a Y Y m LU 
= i — — | جه ص‎ m 
- — - = ( -— - " 4» ** езд» ٤ 


< -4 am —— 


signal depends on the frequency and magnitude of the plant input, R and 
the characteristic of the plant. It is also governed by its original 
signal frequency. If we fix the amplitude and frequency of the input to 
the plant, the attenuation factor is a function of its original perturba- 
tion or drift signal frequency assuming the plant characteristic are con- 
stant. 

From the above discussion, the? is chosen such that the equivalent 
transfer function gives the same phase lag as the plant introduces to the 
same signal coming through the plant input. 

Next, the K is decided by actual measurement of attenuation. 

In measuring the attenuation, the input to the plant is kept constant and 
perturb or drift the parameter by some signal and measure the error out- 
put, which gives the attenuation factor between parameter drift input and 
the error output. Let the measured attenuation be А Тһеп 
^m 7 — Om 
(1 +02 رو‎ 
where Ф n is a parameter drift frequency used 
for measuring А, 
From equation (3-2) K can be calculated. 


The IP curve is represented by IP(X) - ау (К-х) + a, as explained in 


2 


is used as squarer and a, is the coefficient of 


Section 2. A multiplier M 1 


1 


the squarer. After the squarer, а, is added. A proper choice of а, апа а, 
gives an approximate representation of any single dimensional IP curve. 
The bandpass filter has the following form, 


2 S 
E so (3-3) 
S + 2/ 4,5 + (0 


The detector M, is a multiplier. 


2 
9 
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The smoothing filter G, (s) takes various forms depending on the re- 
quirement of compensation of the adaptive loop but generally it is a phase 
lag filter. 

The reference signal filter is an arbitrary filter chosen such that 
phase shift at frequency W due to this filter is equal to the phase 
shift introduced to the perturbation signal between the perturbation in- 
put and detector through the adaptive loop. 

The integrator is usually a D. C. motor and its transfer function is 
designated as G, (s). 

The complete block diagram of a single dimensional parameter perturba- 


tion adaptive control system becomes as Fig. 4. 
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3.2 Modulation 
The input signals through the parameter are the perturbation signal 
and parameter drift signal due to the variation of the parameter value. 


The instant error deviation of parameter, Y (Е) сап Ъе ехргеззе4 аз 


Y, (€) = X(t) +C, - X (3-4) 


O 
where X(t) = instantaneous parameter drift 


Xo = optimum parameter value 


C 
P 


The analysis of signals will be made for two cases of parameter change; 


parameter perturbation signal 


sinusoidal drift of parameter and step drift of parameter. 
(a) Sinusoidal drift of parameter 

The parameter drifts sinvsoidally around the optimum parameter value 
with very low frequency. 51 = C,Sin (0, t is used as perturbation signal 


throughout the analysis. The instantaneous parameter drift is 


X(t) = X, +X, 


sinw t (3-5) 
a = parameter drift frequency and is much lower 
frequency than perturbation frequency, e 


Then Y, (t) becomes as below using equation (3-4), 


Y (t) = Csin Wt + Xi sin (Y t (3-6) 


Referring to the Fig. 4, Y, (t) goes through G (s) and comes out as 
Y,(t), 


1) (O) t - ° E A X sin( (O t - 8) (3-7) 


Y(t) = e(t) = А С 


where A , А are associated attenvation factors, and 


m 9 are associated phase shifts. 


12 





Signals at each point in the loop are labeled in the Fig. 4. 


Following the lcop, 


Y,(t) a, о) 4 (3-8) 


Y, (t) = a, + Y4(t) (3-9) 


Calculating above equations, Y, (t) becomes 


AC AX1 رہ‎ cos2 
Y, (t) = a, + a, ( 5 +0, — - 5 (wt - в) 
222 
7 ا‎ cos2(( t - 0) (3-10) 
A جج‎ m m 


+ sino! - Ө )зіп( 16 - ey 
The last term of Y, (t), a 


modulated signal term where the perturbation signal is amplitude modulated 


sin(w,t - o p sinl wt - 2 is the 


by parameter drift signal. This term carries not only the information of 
magnitude of drift signal but also the information about the actual operat- 
ing point of parameter, in other words, the actual parameter is higher or 
lower than the optimum parameter value by the following form in an actual 
signal. 


When X(t) - Xo O or X sin(@ t - ө) > 0, this term is 
2а AAC, | x,sin(w_t - 0) | sin(w,t - 2 (3-11) 
When X(t) - X, < O or X,sin(4 t - З 0, this term is 


2a.A AC 
m 


1А | Xx sin(w,t - ө) | sin(w,t - ae х) (3:12) 


This is explained in Section 2 and shown in Fig.3. 
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(b) Step drift of parameter 
The parameter value drifts in step from optimum parameter value to 


some value X, then by equation (3-4), the instantaneous deviation of para- 


meter becomes 


Y (t) =X+C - Xo = К-Ж) + С. 511 (3-13) 


Going through the equivalent transfer block, G (s) and comes out as 


Y, (t) selt) = A(X - X.) + жш. - 2 (3-14) 


where Ago 7 are associated attenuation factors and 
58 is the phase shift. 
Calculating the same way as in case (a) using equations (3-8) and (3-9), 


Y, (t) comes out as 


2 2 2 2 
2 23١ SD | 
Y, (t) =a, + 81 Gc - Хо) + 5 J = 5 cosw,t 
+ > - Xp) sin wt - в) (3-15) 


In Y, (t) the last term, — G: - Хо) $11 (42 Е - s an amplitude modulated 
signal term where the carrier is the perturbation sígnal, sin СД: - E 

This term carries the information about the amount of parameter drift (X-X)) 
and furthermore the direction of parameter drift from the optimum parameter 
value. This term takes following form for actual signals: 


When X - Xy ? O or the actual parameter is higher than the optimum 


parameter value, this term becomes 


2а (А АС, | x D | sin((0,t - Ө) (3-16) 


When X - X. < 0 or the actual parameter is lower than Xo» this becomes 


14 





X - X | sin( (2 t - k: - л) (3-17)‏ | اہو سے 


as explained in Section 2 and shown in Fig. 3 
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3.3 Filtering 

The bandpass filter G) (s) attenuates considerably all unnecessary 
signals in Y, (t) and only passes the components that contain information 
about the parameter drift, without much attenuation. 
(a) Taking the case of sinusoidal parameter drift, the У, (=), equation 


(3-10) in Section 3.2(a), is rewritten as 


cos2(w t = om) 


2 


Y, (t) 7 Dg - D,cos2(w,t - ө.) = D 





+ D,sin(w t - 9) 511 (1 Е - БЕ) (3-18) 
where 
В 247 22 
Do а. + > 82 A Xy) 
a Ae 
р. = 1 р 1 
1 2 
а Ағу” 
D. = lm 1 
DR 2 


D, = +7۶ 


Expanding the last term of equation (3-18), 


Y, (t) = 7 - D,cos2(w,t - رت‎ - D,cos2( gy, t - 9 


oO 


D , 
cos 3 сов _ 
[ (ал - Y )t = E: * 3 — [ + (2 )t 8 


+ 
| 
(9 


P 
Е е) (3-19) 
The Y, (t) contains the following frequency signal components, 


| E Е 
Б en ©“ e, — ہت‎ 
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The necessary intelligence is contained in w, = Y and (0, + a com- 


1 
ponents. Therefore the bandpass filter is designed to pass these two 
components without attenuation and suppress other components. Thus the 
center frequency of the bandpass filter should be e and this filter does 
not introduce any phase shift to the signal of ал. Assuming that the 
phase angle characteristic is linear near the center frequency, which is 
generally true for an actual filter, the signal of ( e - رہ‎ gets +A 8 
phase shift and the signal of ( e + w ) gets - A 8 phase shift with 


respect to the signal of W Then Y, (t) passing the filter, becomes as 


1° 


Y; (t) = A Do - A,D,cos2 (Е - 9, - 94) - А, соѕ2( 4) Є - @ - 4) 


A,D, cos 
— 


w - И 
( 1 wW )t с a o) 


Ар; cos 
— С + 0 dt = ° = 8 > ae] (3-20) 
where Ay» Ау, Ay; Аз, and А, are attenuation factors associated eith each 


signal and 


0,, 9, and A O are additional phase shifts introduced to each signal. 


P = 


Usually in a band pass filter A, and А, ате the same attenuation as the 


3 


center frequency signal suffers because WK ал as assumed in Section 3.2. 


Ag» А, апа А, 3 


the last two terms of equation (3-20), Y; (t) becomes 


should be much smaller than A, for our purpose. Combining 


- A.D сов2( 41 - ВЕ - 91) - A,D cos2( t - 8 - 9,) 


75 (Е) = AD - Аур, 20; 


+ A¿D¿sin( wt - o) ns - Ө - A9) (3-21) 
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If the modulated signal term in Y, (t), equation (3-18) is compared with 
the modulated signal in Y (t), equation (3-21), it is moted that the band- 
pass filter introduced a phase lag of Á 0 only to the modulating signal 


component, (sin Wt - Ө) as far as the modulating signal term is concerned. 


(b) In step parameter drift case, Y, (t), equation (3-15) is rewritten as 


Y, (t) р, = 02 w yt + 06 


X - (i - ө.) (3229) 


2و 
h D a, +a A (X -X а + —‏ 
s 2 1 d o) 2‏ 


a Е 
1р1 


6 2 


0 
H 


2а А А С 


Р7 Juri 


Passing the bandpass filter, Y, (t) becomes 


Y.(t) = A.D. - A)Decos(2 w 


106 t - ө.) + A.D 305 (Х - X,)sin(&)t - 7 (3-23) 


1 
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3.4 Demodulation 

The modulated signal У. (Е) reaches the detector Mo; where Ys (t) is 
multiplied by the detection reference signal and demodulation is accomplish- 
ed. The detection reference signal should be in the same phase as the car- 
rier signal to get maximum output at the demodulation. For this purpose, 
detection reference signal phase lag filter G4(s) is used to give the refer- 
ence signal an exactly same phase lag as the phase lag introduced to the 
carrier signal between the perturbation input at the plant and the detector; 
in our case which is ~ Thus the detection reference signal has the follow- 
ing form, 


Y (t) = BC, sin( wt - 05) (3-24) 


where B = attenuation factor of G, (s) 
% = phase lag introduced by 6. (Е) 
(а) For the sinusoidal parameter drift case, Y. (t) equatlon (3-21) is 


multiplied by Ув (Е) equation (3-24) and we get 


У. (Е) = AyD BC, sin( wt - En - А, В C,cos2( Wt - E - 91) 


Е - 22, - A,D BC cos? (w t - em - 9,)sin(w,t - 3 


sin( Ww 2 1 


1 


2 
= i — = =? 
+ 8ھ‎ C,sin (wt AIEN wt e, 40) (3-25) 


In У. (Е) only a D.C. component is integrated by the integrator following and 
develops a parameter correction signal. Only the last term of У. (Е) contains 
a D. C. component which is varying slowing in accordance with the modulat- 


ing signal, which is the parameter drift. All other terms do not contribute 
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to the parameter correction because they are comparatively high frequency 


sinusoids and furthermore, Ао» A, and А) are very small due to the band- 


pass characteristic of G) (s) and they will be further attenuated Ly the 
smoothing filter following. If only the last term of Y, (t) is expanded 


it becomes 


A D4B C 


2 


A D.B C 
373 Isin (:/0,t - 9 - A9) - 


; sin( wt - 8 - AQ) 


cos2 (wW t - > (3-26) 


The first term of (3-26) is considered as a slowly varying D. C. which 


corrects the parameter drift. Thus, 


A,D.B Cy 
You. c. (t) * — 75 —-sin( Om SO A 8) (3-27) 


Putting back the original value for р. and taking into account the phase 


shift between Y. (t) and Yo (t), 


2 ۱ 
Yo (t) = مت‎ Б CiX,sin( Qt - o - А6- ө.) (3-28) 


where 9. is the phase shift in G, (s) and G, (s). 


А, is the attenuation factor associated 
The parameter correction signal Yo (t) has the same form as the parameter 
drift and has some additional phase lag compared to the parameter drift. 
This correction signal will be fed back negatively to parameter adjustment 
and corrects the drifted parameter to the optimum parameter value. 


(b) For the step parameter drift case, Y. (t) equation (3-23) is multiplied 


by Y (t) equation (3-24) and we get 


У. (Е) = A D.B C _ sin( w t - => - A,D,B C cos(2W,t - 9,) 


0 5 1 16 1 
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T2 
+ A¿D,B с, (Хх - Ko) sin ( wt - > (3=29) 


In У. (Е) only the last term develops a D. C. component, which will be 

integrated and drive the parameter adjustment. Other terms do not contri- 
bute to parameter correction because they sre comparatively high frequency 
sinusoids and Ay and Ay are very small due to the bandpass filter attenua- 


tion and will be further attenvated by the smoothing filter. If only the 


last term of Y, (t) equation (3-29) is expanded it becomes as, 


1 
A_D_B C, (X-X_) cos2( ty, t - 98) (3-30) 


"0 ۵0 


2 3% 
Putting back the original value for D,, 


2 
Ya, c, C? = а.А mo T ED (3-31) 


ld o? 
Thus the rate at which the parameter drift is corrected, is proportional to 
the parameter error (X-X9), the perturbation signal ee С,» the con- 
stant а, associated with the IP curve (the sensitivity of the IP variation 
to changes in parameter value) and adaptive loop gain. This open loop 
analysis is not exactly correct when the loop is closed, but it shows a 
possible signal development through the acaptive loop. 
(c) Next, the effect of phase difference between the carrier and the detec- 
tion reference signal Y (t), in the demodulacion, is considered, To see 
more clearly the process of demodulation and the effect of the phase dif- 
ference between carrier and the detection reference signal, graphical multi- 
plication of actual signals involved, is performed, Taking the case of 


step parameter drift, the actual modulated signal in Y. (t) has the following 


forms as explained in 3.2(b) and expressed in equations (3-16) and (3-17), 
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а |х - xo | sin(w,t - 8) when x- x > 09 (3-32) 


A| x - x | sin(w,t - 8, - m) when X- X0 (3995) 


where À is a constant, 


The carrier in У. (Е) is C(t) = sin( Qt - РВ) when X-X4 2 0 and -€, 5 (t) = 


sin(w,t - 98 - x) when Х- Хо < 0. 


Тһе detection reference signal is Y (t) = B C sin( Wt - 9). 


1 
If e. Ё SS Y (t) leads or lags the carrier E by some angle. The 
graphical multiplication is performed for the following three cases, 

(1) Y (t) is in phase with the carrier. 

(2) Ув (Е) leads the carrier by 30°. 

(3) Y (t) lags the carrier by 30°. 
The coefficients, of each term is a constant, therefore only a sine wave of 
unit maximum amplitude represents each wave. 

Fig. 5 shows the demodulated output by dotted line for the above three 
cases. It shows clearly the phase MIU nature of perturbation signal; 
when the carrier in a modulated signal is EU = біп (42,6 - nr the 
demodulated output is all positive and when the carrier is Cs = sin 
( wt - E - X) the demodulated output is ail negative. Thus it will give 
correct direction of parameter correction in accordance with the actual 
parameter error, X 2 О ок X-X4 L 0. 

In case (1), we get maximum net output of correction signal when inte- 
grated as seen in Fig. 5(1). For the case (2), we get less net output of 
correction signal when integrated. This eventually ends up with less adap- 
tive loop gain. The case (3) has the same output as in the case (2). 


If we consider the time lag of starting to build up correction signal 
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output at the integrator, the case (3) has the largest time lag of the 
three cases, and the case (1) is the second and the case (2) is the least. 
But this time lag is negligible because the carrier frequency is very high 
0 to the modulating frequency. 

From this analysis, we can see that when the phase difference between 
Ys and Ye is exactly 90°, the signal output will be the lowest and when the 
phase difference is more than 90° the demodulated signal polarity will be 
reversed and the magnitude will increase as the phase difference approaches 
180°. Therefore, when the carrier phase lag is expected to be closer to 180°, 
it is better to remove the reference signal phase lag filter and change the 


polarity of the demodulated output to have correct feedback. 
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4. The Simulation Study of a Sinusoidal Perturbation System 
4.1 IP. variation investigation 

To investigate the actual variation of the IP. with respect to the 
parameter drift and perturbation of the plant, the following model and 
plant were simulated by an analog computer with the adaptive loop opened. 
The analcg computer used throught the investigation is the Donner 3100. 


K 


Model; Gnog lE? = ECE ES 


K 


Plant; G (s) = S(S^ x) 1 


The pole of the plant is a variable parameter x, and the pole of the model 
is fixed at 3, which is the optimum parameter value in this case and gives 
Ej of 0.5. The pole of the plant can vary in the range that gives of of 
plant from 0.1 to 1. with K= 9. The block diagram of the simulation is as 
Fig. 6 and the analog computer circuit is shown in Fig. 7. 

(a) The input to the plant and model is a square wave of 0.1 cps. 
A plant parameter, the pole of the plant is drifted by a triangular wave 
of 0.01 cps. while the model and plant are operating. The output of model 
and plant are as Fig. 8(a) and (b) respectively. The output wave shape 
shows the change of f of the plant due to the parameter drift. 

(b) The error output is recorded as Fig. 8(c) and the pole drift 
is as Fig. 8(d). 

Next, the input to plant and model is sine wave of C.3 cps. and the 
error output is also recorded as Fig. (а). The Fig. 9(b) is the pole drift 
as before and (c) is the error squared output. 

The above results shows that error is produced only when the output 
of plant is in transient and the plant response is nct the same as the 


model response due to the parameter drift, Also it shows that the magnitude 
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Fig. 8. The output of model, plant and error.with plant pole drifted. 
a) model, 


(c) error 
8 plant pole drift 
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Pig, 9. Error and error squared output with sine wave input to plant. 

a) error 
R plant pole drift 
(c) error squared 
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(c) 





of error is proportional to the parameter error, the magnitude of para- 
meter drift. 

The magnitude of error is bigger when the parameter drift is in the 
direction of decreasing plant J with the same amount of parameter drift. 
Therefore, it is expected that the magnitude of parameter correction signal 
will be greater for the parameter drift to a decreasing p than a increasing 
F , With the same amount of drift. 

(c) Error output is observed with only perturbing the pole of 
plant by a sinusoidal perturbation signal. The perturbation frequency is 
much higher than a drift frequency; here a sine wave of 5 cps is used, 

(1) The magnitude of error output due to perturbation signal 
only. The perturbation wave appears in the error output, but the magnitude 
is much reduced due to the severe attenuation, Five cps is a high frequency 
for this plant. Thus the parameter perturbation signal does not appear con- 
spicuously in the plant output even with a big amplitude perturbation signal 
as far as the frequency of perturbation is high. ТЕ we want to eliminate 
the effect of perturbation on plant output completely, the parameter of 
model can be perturbed instead of plant parameter. Fig. 10(a) is the per- 
turbation wave appearing on the output of the plant. The perturbation 
waves appear conspicuously on the output of the plant only when the ratio of 
the amplitude of the plant input to the amplitude of the perturbation signa! 
is small, 

(2) Phase shift of perturbation signal through plant. 

The perturbation wave appearing on the output of plent lags the original 
perturbation signal by some angle. This phase lag is introduced by the 
plant and depends on the characteristic of plant and signal frequency. 


Fig. 10(c) and (d) are the original perturbation signal and the perturbatio: 
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(c) Original perturbation. P Input to plant 
(d) perturbation wave on output (f) Output of plant 


Fig. 10. Perturbation waves on the output of plant. 
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wave appearing on the plant output respectively. The same marks on the 
waves are the corresponding points on both waves. These corresponding 
points are found by starting and stopping the signal abruptly. According 

to the measurement of phase difference between the two waves, perturbation 
wave appearing on the output lags the original signal by 165°. To compare 
this phase lag with the phase lag introduced by plant to a signal coming 
into the plant through the plant input, the same signal used as perturbation 
is fed into the plant. Fig. 10(e) and (f) are the input to plant and the 
output of plant. The same marks on the signals are the corresponding points 
on both waves which is found as before. According to the measurement, the 
output waves lags input by 165° approximately, which is the same phase lag 
introduced by plant to the parameter perturbation signal. From this re- 
sult it is concluded that the phase lag introduced by the plant to the 
parameter perturbation signal fed through the pole of plant is the same 
phase lag introduced by plant to a signal of same frequency fed through the 


plant input. 
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4.2 The Simplest Adaptive System Simulation (Adaptive circuit 1) 

The complete adaptive control system is as Fig. 1 and the block dia- 
gram is as Fig. 4. To investigate the function of each component and the 
effect of change in components on the over-all system operation, the simula- 
tion is started with the simplest combination of components for a adaptive 
system. Then more components are added to improve the system operation 
until the complete adaptive system is formed. We designate successive cir- 
cuits by number from the simplest to the complete circuit in simulation 
study. The first circuit is designated as adaptive circuit 1, which is the 
simplest and is shown in Fig. 11(a). The analog computer set-up is as Fig. 
11 (b). It simulates starting from the error output point in the complete 
system, In Section 3.2, signal analysis, the instant deviation of parameter 


Y (t) is expressed as equation (3-4), 


Y, (t) = X(t) + са - X (4-1) 


0 


For the analysis the optimum parameter value X, is arbitrary and can be 


0 
set as zero. Y (t) causes an error output and the error output Y, (t) has 
severe attenuation and phase shift and comes out as equation (3-6) for a 


sinusoidal parameter drift. But in the simulation study, ignore the at- 


tenuation and phase shift and assume the error output Y, (t) as, 
Y, (t) = X(t) + E (4-2) 


This error output is going to be the input in the Adaptive circuit 1 because 
the circuit simulates from error output point. In an actual system, the 
parameter correction signal corrects plant parameter and reduces the error 
output to zero but in this simulation the correction signal directly reduces 


the У. (Е), which is equivalent to correcting the plant parameter. 
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Fig. ll. Adaptive circuit 1. 





In the circuit, м апа M, are a squarer and a detector, both of which 
are electronic multipliers. Тһе zero D. C. level of electronic multiplier 
is not constant, therefore a pure integrator can not be put in right after 


the multiplier M, because the integrator integrates the drifting D. C. and 


2 
K 
developes erratic correction signal. A pseudo integrator — — is used 
K ۰ 
for the simulation of the pure integrator — . The parameter perturba- 


tion and parameter drift signals used are as follows: 


5 cps. 


сЕ = C‚sinw,t 2 ЕТТ (2, 


КСЕ) =X sinw 6 Х. = 8+. Ш - Oclwcpa. 


1 т 


feedback point and open loop observation is 


1 
The adaptive loop is open at Y 


9 
made, 

(a) The wave shapes of signals at every point of loop are shown 
in Fig. 12. The (a) in the Fig. is the error signal in which perturbation 
signal is superimposed on the parameter drift signal. (b) is the output of 
squarer where the modulation is accomplished. (c) is the output of detector; 
the demodulated output has positive or negative D.C. component in accordance 
with the parameter drift direction. This is integrated and comes out as (d), 
which is the parameter correction signal. If the parameter correction sig- 
nal is compared with the parameter drift signal in (a), it has the same shape 
as the parameter drift, therefore, when the correction signal is negatively 
fed back it will correct the parameter drift. But the correction signal lags 
the drift signal by 60? approximately as seen in (a) and (d). The marked 
points are the corresponding point on both waves. This phase lag is caused 
by the integrator as seen between (c) and (d). There is no phase lag be- 
tween Y, and Y, as seen in (a), (b) and (c). The perturbation signal ap- 


1 7 


pears in the correction signal Y,, which will be eliminated by adding a 


9? 
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smoothing filter. 
A square wave also can be used as a parameter perturbation и 
The same results as obtained by sinusoidal perturbation, are obtained 
using square wave perturbation. The square wave perturbation wave ap- 
pearing in the correction signal looks like a sinusoidal wave. Fig. 12-1 
shows the results. 
(b) When the frequency of perturbation is increased, it is at- 
tenuated more and it does not appear in the correction signal; in this 
case there is no perturbation signal in the correction signal when 15 cps. 
is used for 421. The amplitude of the correction signal is not affected 
by the change of perturbation signal frequency. This means the perturba- 
tion frequency does not affect the loop gain in this simple circuit, how- 
ever, in an actual circuit the frequency of perturbation will affect the 
amount of attenuation due to circuit components and affect the loop gain. 
Fig. 13 shows the correction signal with different perturbation frequencies. 
(c) Next, the perturbation frequency is fixed at 5 cps. and the 
frequency of parameter drift is varied. When the drift frequency increases, 
the amplitude of correction signal is decreased, Fig. 14 is the amplitude 
of correction signal vs. drift frequency graph. The decrease of amplitude 
of correction signal is due to the attenuation by the integrator. This re- 
sult and the result from (b) show that in a closed loop analysis of the 
adaptive circuit, we need to deal with the parameter drift signal only. 
The perturbation signal does not affect the magnitude of correction signal 


even though it comes through the detector and rides on the correction signal. 
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Fig. 13. Correction signal with different perturbation frequency. 
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4.3 Adaptive circuit 2 

From the simulation study of adaptive circuit 1, it is found that a 
smoothing filter is needed to eliminate perturbation signal appearing in 
correction signal. A smoothing filter G (s) is inserted between M, and 
the integrator in adaptive circuit 1. The circuit becomes as Fig. 15(a). 
The analog computer circuit is as Fig. 15(b). This circuit is designated 
as adaptive circuit 2, The smoothing filter will filter out completely the 
perturbation signal appearing on the correction signel but also attenuates 
the correction signal and gives more phase lag to the correction signal, 


Too much phase lag of the correction signal tends to mexe a closed loop 


e 


system unstable, It causes time laz of the correcting parameter and when 
the parameter change is abrupt such as step change, excessive overshoot re- 
sults on the correction signal and sometimes it produces a sharp undershcot 
before producing correction signal in correct direction, The undershoot is 
a strong correction signal in à wrong direction. It will be studied in more 


detail later. The smoothing filter is usually chosen as a compensator for 


the adaptive loop to improve step response of the system. 





(a) In the simulation, a filter of the om, T3 S is used. 
The E values used are 0.2, 0.5, 1, 2, 5, 10. The same perturbation signal 
and — drift signal used in 4.2, are used, Taking the perturbation 
signal present in correction signal as noise, correction signal to noise ratio 
is measured at the output point of the smoothing filter Co (83e Fig. 15 is 


the graph of correction signal to noise ratio vs. "b, of the smoothing filter. 
< 


According to the graph, Ta of 2 or above gives more than 10 : 1 signal to 
noise ratio. The perturbation waves remaining in the correction signal even 


after the filtering, will be further attenuated by the integrator following. 


c of 2 gives corner frequency of 0.5 rad./sec, erá this gives enough 
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(b) Computer circuit 


Fig. 15. Adaptive Circuit 2. 
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attenuation to perturbation waves of 5 cps. present in correction signal of 
this system. Thus it is concluded that an acceptable ratio of corner fre- 
quency of smoothing filter to perturbation frequency is 1 : 60 or greater. 
But the attenuation of perturbation waves in the correction signal is not so 
important as the stability of the system, therefore, this filter is designed 
as an appropriate compensator for the adaptive loop and it naturally acts 
as a smoothing filter, 

(b) Effect of adaptive loop gain on parameter drift correction 

(1) Sinusoidal parameter drift 
The parameter is drifted sinusoidally. When the adap- 
tive loop is open the parameter drifts sinusoidally and as soon as the 
adaptive loop is connected the drifting parameter is corrected to the optimum 
value. The loop gain affects the amount of steady state error as expected. 
The necessary loop gain is decided by the required steady stete accuracy. 
Fig. 17(a) shows a parameter drift when the adaptive loop is open and then 
corrected parameter condition when the loop is closed. (b) is the corres- 
ponding parameter correction signal. The loop gain is a little low and the 
steady state error is comparatively large. To reduce the steady state error, 
loop gain is increased and the parameter drift is almost completely correct- 
ed as seen in Fig. 17(c). With this high loop gain, the loop was closed 
suddenly when the parameter drift was a maximum, the parameter drift was 
corrected to the optimum value with an excessive oscillation before settling 
down to steady state; Fig. 18(a) shows the parameter drift being corrected 
and (b) is the corresponding correction signal wave shape. 
(2) Triangular Parameter Drift 
Fig. 18(c) and (d) is the triangular parameter drift being 


corrected and its correction signal. The loop gain is a little low and the 
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Fig. 17. Parameter drift correction 
a) Parameter drift being corrected with low loop gain 
b) Parameter correction signal 
(c) Parameter drift being corrected with a' high’ loop:igain 
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Fig. 18. Parameter drift correction 
(a) Sinusoidal parameter drift being corrected with high gain 
(b) Correction signal 
Я o Triangular parameter drift being corrected 

d) Correction signal 
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correction is not perfect. 
(3) Square wave parameter drift 
The parameter is drifted by a square wave. Fig. 19(a) 

shows the square wave parameter drift being corrected and (b) is the corres- 
ponding correcting signal, in which steady state error is small. The system 
is stable. When the loop gain is too high, the system response becomes 
oscillatory and tends to be unstable but the system will never be unstable 
because overall loop is second order system in this set-up. Fig. 19(c) 
is the parameter being corrected and (d) is the corresponding correction 
signal when the loop gain is very high. 

(c) The effect of perturbation frequency on adaptive loop step 
response. 

The most extreme case of parameter drift is step change of 
parameter due to an abrupt environmental change. Therefore, only the effect 
of perturbation frequency on adaptive loop step response is investigated. 

The plant parameter is drifted by square wave of 0.01 cps. For dif- 
ferent frequencies of perturbation, maximum oveishoot and rise time are 
measured, The perturbation frequencies are varied from 1 cps. to 100 cps. 

(1) Steady state magnitude of correction signal 

In this adaptive circuit there is no attenuating com- 
ponent for the perturbation signal between perturbation signal input and the 
demodulation, therefore, the perturbation frequency does not affect the loop 
gain and the steady state magnitude of correction signal is constant with 
different perturbation frequency as seen in Fig. 13. 

(2) Maximum overshoot and undershoot 

For a very low frequency of perturbation, the correction 


signal has extremely irregular wave shapes of square waves as shown in Fig. 


47 





— 


0 ЕГЕЕЕЕЕЕЕЕЕЕЕЕТЕЕЕЕЕЕЕЕЕЕ 


MH EEE EEE 
Pig و کم لد‎ il Г 
аа а. 


-e mi шшш cn 
„Eee کی‎ 

ас у ск 

=з оа э. 


ORDER MARK. 
Е s s 
EEEEEEEEEEEEEEEEFEEEEEEEE = 


“л 





о 


— 





— ہے‎ — — EE 
o EEEE 200 77 
5 Е ІП ЖЕСЕ 
(с) 
—— yep qp Hr ШІП) 
— Sal 


| = = 
0د‎ Eme Eee ae 
= EE EEE کات تحت ا‎ 


— m —— 
20 —— — — ا‎ = 





















во Еш ее == е == 
— е Е Е 
ЕЕЕ 
Sl мацы да 


Fig. 19. Parameter drift correction 
(a) Square wave parameter drift being corrected 
(b) Corresponding correction signal 
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Fig. 20 Parameter drift correction with low frequency of perturbation 
(perturbation frequency is 2 cps,) 


(a) parameter dríft being corrected. 
corresponding correction signal. 
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20(b). The perturbation frequency is 2 cps. The correction signal has 
various amounts of overshoots and undershoots. Here, undershoot is de- 
fined as a sharp and strong correction signal to a wrong direction. The 
magnitude of overshoot and undershoot depends on the point of perturba- 
tion signal wave at which an abrupt change of parameter occurs. The ir- 
regularities of correction signal wave shapes becomes greater at lower 
perturbation frequency and with perturbation frequency below 1 cps. the 
correction signal is too irregular to have a sensible correction, The 
magnitude of overshoot and undershoot also decreases as the perturbation 
frequency increases and above 5 cps. of perturbation frequency there is 
no overshoot nor undershoot. Table 4.1 lists the measured ratio of maxi- 
mum overshoot and undershoot to steady state correction signal amplitude at 


lower perturbation frequencies. 


4 (cps.) Max. overshoot Max. undershoot 
1 1.63 4.25 
1.5 1.28 2.06 
2 1.19 1.81 
3 no overshoot 1.38 
4 ч 1.18 
5 1.03 





Table 4.1 Ratio of Maximum overshoot and undershoot to the amplitude 
of steady state correction signal. 


It is noted from Table 4.1 that the maximum undershoot ratio is much great- 
er than overshoot ratio. The sharp undershoot should be avoided, which 
makes the system unstable and leads to erratic operation. Therefore, the 
perturbation frequency should be higher than 5 cps. in this adaptive control 
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(3) Rise time 
Rise time vs. perturbation frequency is plotted as Fig. 
21. The rise time increases sharply from e = 1 срз. to а, = 4 cps. 
Above 5 cps, rise time increases very little and Stays almost constant. 
Fast rise is desirable and it is obtained at lower a but it accompanies 
undesirable undershoot as seen in Table 4.1. In the following simulation 
studies, e = 5 cps. is used as perturbation frequency throughout, which is 


the lowest perturbation frequency that does not accompany any overshoot or 


undershoot in correction signal. 
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4.4 Adaptive Circuit 3 

A band pass filter G) (s) was added between the squarer and the de- 
tector in adaptive circuit 2, Fig. 15. The circuit becomes as shown in 
Fig. 22. The band pass filter attenuates all unnecessary frequency compon- 
ents considerably in the signal coming from squarer and passes the components 
which contain information without attenuation. The band pass filter is es- 
sential in a multi-loop adaptive control system where different perturbation 
frequencies carry information for different parameters being controlled through 
different loops. This filter also affects the stability of the adaptive loop. 
The requirement of the band pass filter is explained in Section 3.3. To meet 
the requirement a filter of following form is used, 


obs 


(4-3) 


Сс (у= с L 

1 2 2 

S + 2f Ws +4, 
where (4 = perturbation frequency (center frequency) 
2) = damping ratio 
The half power frequency of this filter is calculated as, 

U = ро ЕА 

= 201 4 )4-4( 


The band width %b.w. = 2 f 01 
To decide the f value which gives best step response to the system, the Y 
is varied and the step response of the system is observed. Wi of 0.5 or 
0.6 gives best response, 

(a) First, only the band pass filter is simulated. f of 0.5 is taken 


and the filter becomes as, 





G (s) = І (4-5) 














2 
wS 


5° +2] 05+ а? 


с, (в) с. (в) 
Fig. 22. Adaptive circuit 3 


Y3 
Ys 


Fig. 23. Analog computer simulation of band pass filter с. (в) 
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The analog computer simulation of the band pass filter is shown in Fig. 23. 
The frequency response and phase angle characteristic is plotted as Fig. 24. 
The calculated half-power frequency and band width and the actual half- 


power frequency and band width are as below, 





(A, (rad/sec) ал „(таа/ѕес) Wy, wy, (rad/sec) 
calculated 15.7 &7.1 31.4 
actual 19.2 49.5 SO 





Using the result of signal analysis in Section 3.3(a), following signal 
components are expected to be present in the signal from the squarer of 


adaptive circuit 3. 


2С 
20 = 1,256 rad/sec 
244 - 62,8 р 
e - w = 30.144 rad/sec 
— " 
ал % «АЗ 32.656 


Therefore, this filter will suppresses D.C., 20, апа 2 e, and passes 


CU, - i, and w+ 0 a components without attenuation. 


(b) The adaptive circuit 3 simulation was obtained by connecting 
the band pass filter simulation, Fig. 24 between M and M, in adaptive 
circuit 2, Fig. 15. The performance of the adaptive loop is much improved 
by the addition of the band pass filter. With proper loop gain, a square 
wave parameter drift is corrected perfectly as shown in Fig. 25, where (a) 
is the square wave parameter drift being corrected and (b) is its correction 
signal. The ШЕ is 1.08 and the rise time is 0.7 sec. Even with excessive 
loop gain, the performance is less oscillatory compared to the same case in 


adaptive circuit 2, in 4.3(b)(3). Fig. 25(c) is the parameter drift being 
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Fig. 25. Parameter drift correction 
(a) Square wave parameter drift being corrected 
b) Corresponding correction signal 
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corrected with excessive loop gain and (d) is its correction signal. 
Compare these figures with Fig. 19(c) and (d) which are the same case in 
adaptive circuit 2, Fig. 25(c) shows that the parameter drift is over- 
corrected because of excessive loop gain but the response is less oscil- 
latory. The comparison shows a definite improvement in adaptive loop per- 
formance by adding band pass filter. This filter is acting as a compensator 
in this circuit. 
(c) A comparison of demodulated outputs is made between following two 
cases in Sinusoidal parameter drift. 
(1) Modulated wave is not filtered by band pass filter. (adaptive 
circuit 2). 
(2) Modulated wave is filtered by band pass filter. (adaptive 
circuit 3). 
Fig. 26(a) is for the case l. Refer to the Fig. 22 adaptive circuit 3 
for the designation of signals at each point. For the case 1 the modulated 


signal Y, goes into the detector without filtering by band pass filter. 


3 


The demodulated output is Y The modulated signal Y, contains many un- 


7° 3 
necessary signal components which do not contain any intelligence as ex- 
plained in section 3.2. Accordingly, the demodulated output contains un- 
necessary signal components., On the other hand, in the case 2, Fig. 26(b), 
the У. passes the band pass filter which filters out unnecessary signal 


components as explained in section 3.3. Thus, the signal Үс contains only 
signals which contain intelligence. Note the wave shapes after filtering. 
This signal is demodulated and comes out as Y, which has clear.cut polarity 
in accordance with the parameter drift direction. This indicates the Y, 


for the case 2 does not contain any unnecessary signal as in the case l. 


This verifies the signal analysis in section 3.3 and 3.4. 
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(а) 


and unfiltered, 
(a) Unfiltered 








4.5 Adaptive circuit 4 

The measurement lag filter G(s) is connected between the parameter 
perturbation and drift input and squarer. The detection reference phase 
lag filter G4(s) is connected between reference signal input and the de- 
tector. The G(s) represents the phase lag in measuring error. 6. (3) 1$ 
chosen to give the same phase lag to the detection reference signal аса) 


as G(s) gives to the perturbation signal (the carrier) at &,, thus the | 


1 3 


carrier and the detection reference signal are in phase with each other at 
the detector. The band pass filter G. (s) does not introduce any phase lag 


tn the carrier at eo, and the total phase lag introduced to carrier be- 


tween Y, and Y,, is the only phase lag due to G. (s). The circuit becomes 


1 5° 


as Fig. 27 and analog computer simulation is as Fig. 28. This circuit is 
designated as adaptive circuit 4. The filters G(s) and 6. (3) take follow- 


ing forms, 


ne т т i G(s) 1+ 2.5 


(а) The process of demodulation is explained and analyzed in section 


1 
. e e е е е е е am * Sa eee 
3.4. It is verified in this simulation study. First, ¿= Ç= 31.4 


7, 


is used, which gives 45° phase lag at “| - 31.4 rad/sec/(Scps.). Thus 


the carrier signal in Y. and the detection reference signal Y, are in phase 


6 
is multiplied by the Y 


5 


at the detector, The modulated signal Y signal 


5 6 


and the demodulated output is produced. The demodulated output has one 
polarity or the other in accordance with the parameter drift direction with 
respect to the optimum parameter value. Fig. 29(a) is the modulated signal 


Ү. апа (b) is the demodulated output Y The demodulated output has exactly 


5 7 


the same wave shape as obtained by graphical multiplication of Үс and Ye in 
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section 3.4 Fig. 5. This signal is integrated and develops a correction 
signal of one polarity or the other. When the polarity of detection re- 
ference signal Y is changed, the demodulated output also changes its 
polarity. These results confirm the explanation of the demodulation in 
section 3,4. 

The modulated signal and demodulated output wave shapes for square 
wave parameter drift are shown in Fig. 30. The (a) in this figure is para- 
meter drift, (b) is the output of squarer, (c) is the modulated signal com- 
ing out of the band pass filter and (d) is the demodulated output. The 
demodulated output has one polarity or the other depending on the direction 
of parameter drift. 

(b) Next, the effect of phase difference between Ys and Y, on the 
demodulated output is observed. This effect is explained in section 3.4 
and the result of graphical multiplication of two signals which are out of 
phase with each other by some angle is shown in Fig. 5. The response of 
the system is observed for the following three cases. 

1. Y, and Y, are in phase. 


6 5 


2: Ye leads Үс 


3. Ye lags Ys 


For the case one, the parameter drift being corrected and parameter cor- 


by 60°. 


by 60°. 


rection signal wave shapes are as Fig. 3l(a) and (b) respectively; the one 
is 1.046 and rise time is 0.9 sec. For the case 2 and 3, the results are 
the same and the result of the case 2 is shown in Fig. 3l(c) and (d). The 
response is overdamped and the rise time is about 2 sec. This shows that 
the overall loop gain is reduced due to the phase difference between the 
Y. and Y,. This effect is predicted in the signal analysis in section 3.4. 


5 6 


The phase difference between Ys and Ye does not affect the system seriously 


and merely reduces the loop gain which can be easily restored. 
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Fig. 30. Modulation and demodulation of square wave parameter drift. 
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Fig. 31. Effect of phase difference between Y. and Y. 
(a), (b) Y¿ end Y, are in phase. 
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4.6 Complete adaptive control circuit (adaptive circuit 5) 

Finally, a complete adaptive control circuit was simulated on an 
analog computer. The block diagram of the complete system is as Fig. 32. 
This circuit will be designated as adaptive circuit 5. The analog computer 
simulation is achieved by merely replacing the G (8) in the adaptive circuit 
4 simulation, Fig. 28 with the model and plant simulation circuit in Fig. 7 
and the correction signal is negatively fed back to the summer, output of 


which goes to the multiplier M The reference signal phase lag filter 


3° 
G,(s) is not used because the carrier is expected to have phase shift of 
160° at the detector, which is nearer to 180°. The input to the plant and 
model, is R(t) = 7sinl.884t (f = 0.3 cps.). 

(a) The smoothing filter G, (s) is acting aS a compensator. To see the 
effect of the G, (s) on the system response, first it was removed and the 
response of the system was observed. The adaptive loop was opened at the 
output of the integrator G, (s). The parameter correction signal developed 
with the loop open, is recorded as Fig. 33(b). Fig. 33(a) is the square 
wave parameter drift signal. To close the adaptive loop, the polarity of 
the correction signal should be reversed to have negative feedback, 

When the loop is closed, this system becomes unstable, so the loop gain 
was reduced to a very small value such that the system barely maintains a 
stable condition. The correction signal becomes as shown in Fig. 33(c). 

The magnitude of the correction signal becomes very small and the signal 
shape is drastically changed compared to the correction signal shape with 
the loop open. That is why the open loop signal analysis is not exactly 
valid when the loop is closed as mentioned in section 3. Fig. 33(d) is 
the plant output and (e) is the model output. Because of too small loop 


gain and accordingly too small correction signal, a slight correction of 
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parameter drift is achieved. Compare the figures (d) and (e). When the 
plant parameter changes abruptly the plant output has sharp overshoot due to 
instability of the system. 
1 
(b) A lag filter of the form G, (s) EDEN TS 
tween the detector and the integrator. Three different values of 75, 2, 


is inserted be- 


5, and 10 are used and the loop gain is adjusted to give stable operation. 
With different values of ¢,, the correction signal wave shapes are differ- 
ent but in all cases the critical loop gain is too low. Thus the para- 
meter drift can not be corrected completely in this circuit. Fig. 34 is 
the parameter drift being corrected and corresponding correction signal 

for the three cases; Fig 34(a) is for 7, =12 (b) is for T, = 10, (c) is 
for CT, = 5. An effective compensator must be found to achieve a re- 


quired performance of the system. 
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5. Design of a Sinusoidal Perturbation Adaptive Control System 

A design of a sinusoidal perturbation system is attempted based on the 
simulation study and using a linearization technique. The design procedure 
and considerations are virtually the same for this system as general control 
systems except some special features such as modulation and demodulation of 
error signal using perturbation signal. Finding the approximate equivalent 
transfer function of circuit components and applying linearization technique, 
a preliminary design is made on analytical basis and the preliminary design 
is simulated on analog computer and the effectiveness of the analytical 
method is checked. The final refinement of design should be done on experi- 
mental basis. 
5.1 Specification 

There are many ways to specify the performance criteria. Those speci- 
fications differ depending on the application, but at least the steady state 
requirement and dynamic or transient requirement should be specified. One 
special specification for this system is hunting loss, wnich is defined as 
the maximum allowable deviation of parameter caused by sinusoidal perturba- 
tion at the optimum parameter value. The choice of the hunting loss depends 
on the requirement of the system and the characteristic of the IP curve; a 
small hunting loss gives less distrubance on the output of the plant but may 
make the system lock into small local minimum in IP curve and sometimes 
does not give any information about the actual parameter drift because of 
attenuation. 
5.2 Selection of Perturbation frequency 

The perturbation frequency ал should be much higher than the expected 
maximum frequency of the parameter drift АҒ this is necessary in any 


modulation system. The investigation of the effect of perturbation frequency 
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on the system response in 4.3(b) shows that a perturbation frequency of 
higher than 5 cps, does not make much difference on the system response. 
At an extremely low frequency of perturbation, the system operation tends 
to become erratic due to sharp overshoot and undershoot even though the 
rise time is very short. If the perturbation frequency is too high, the 
attenuation through the system components is severe and the higher gain of 
perturbation is required. 
5.3 Linearization of the System 

To investigate the stability and response of the system analytically, 
an equivalent signal transfer function of the adaptive loop is established. 
The representative circuit block diagram is derived in section 3.1 and 
shown in Fig. 4. The following linearization technique is applied to 
linearize the adaptive loop. In the adaptive loop, a D.C. parameter correc- 
tion signal is proportional to the plant parameter error X - Хо» as shown in 
equation (3-31) section 3.4(b). This indicates a possibility of derivation 
of an equivalent linear loop. In applying the linearizing technique to any 
specific system, only the variation of the IP curve about the optimum is 
important. Therefore in the general IP curve equation, set Хо = 0 апа ay = 
0. Referring to the Fig. 27 adaptive circuit 4, which is a complete block 
diagram of the system, the portion between error output and the تل‎ 18 
a carrier system. To linearize the squarer and the detector, the describirg 
function or a first term of Taylor series approximation is used and to get 
an equivalent signal transfer function for band pass filter G) (s) frequency 
transformation technique is used. (Ref. 1) 

The perturbation signal does not get into the p:cture when a closed loop 


analysis is made as verified in simulation study 4.2. Therefore, only the 
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modulating signal (parameter drift signal) is acting in the closed loop. 
(a) The parameter drift signal acts directly on the G(s) and no 
linearization is necessary for G. (s). 
(b) Linearization of M 


1 


Standard describing function technique is used. The input to м, is the 


error output У. (Е) and can be rewritten from equation (3-7) as, 
Y, (t) = oa Wt - ° + A X ,sin( 2 t - 9 (5-1) 


The output from М, 15 Y, (t) equation (3-10). In this equation only the 
last term is the carrier modulated by the intelligence signal and all 
other terms does not contain any intelligence and will be suppressed by 
the band pass filter. Therefore, as far as the modulating signal is 
concerned, only the last term of the Y, (t) is the output from M. The 


last term of the Y, (Е) is written from equation (3-10), 


2ajA,A.C,X,sin( Wt - 6 )sin (wt - 8.) (5-2) 


1 р 1 
The equivalent multiplier gain is the ratio of the peak amplitudes of 
signal input and output terms. Thus, the peak amplitude of (5-2) is 
divided by the peak amplitudeof the modulating signal term in (5-1) and 


one obtains 


Gul = 22007 (5-3) 


This same gain also can be obtained by a Taylor series approximation. 
(c) Linearized equivalent transfer function of band pass filter G, (s). 
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The steady state phase shift due to G) (s) at S = j ,صا‎ is O° 
The linearization is based on the following assumptions; 


1. Phase shift curve is linear at S = j W: 


2. Attenuation at S = 304, + CO) is almost the same as at S - j«,. 
Considering only the modulating signal component, the input to the band 
pass filter is the signal (5-2) and the output is the last term signal in 


У. (Е) equation (3-21). The output is written as, 


AS ЗЕ (Qt - 8,11 wt ~ Ө - А Ө) (5-4) 


The input (5-2) is compared with the output (5-4), the modulating signal 


is attenuated by the factor of A, and phase lag of AO is introduced by 


3 


the G, (s). An equivalent transfer function may now be defined such as to 


approximate the phase characteristics of the original network or, 


A 


=— — —— 5- 
1 % 74 — (5-5) 


© E 


In this case, A is A, which is the attenuation factor of the band pass 


3 


filter асал. and Ti eq is determined by choosing a representative 


value of and С is calculated such that G 


ed 1۴ (s) gives the same 


phase lag at (0 as 40, which results in the original system. 
m 
(d) Equivalent transfer function for M. 


The detection reference signal input to М, 1$ 


Y (t) = 80681٤ - ° (5-6) 
and the modulated signal input to M, is the signal term (5-4). 


The output of M, is the Ya a (t) equation (3-27) and is rewritten as, 


1 
У7а.с. (=) == (2a, А, А E 1 A,)BC, sin( w — о. - A9) (5-7) 


The ratio of the peak amplitude of input (5-4) to output (5-7) for м, 18 
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(5-8) 


Gy2 1 


> BC 
The components after the M, in the adaptive loop act upon the signal terms 
directly, so no further transforms are required. 
5.4 Stability investigation of uncompensated system and compensation 

After linearizing the required شر‎ е. the block diagram of the 
system is formed with only the fixed compcnents. All components are fixed 
except the smoothing filter G, (s) which is acting also as a compensator. 
The block diagram is reduced to a convenient form and the overall loop 
transfer function is obtained, The gain requirement is decided on the basis 
of static and dynamic accuracy and the gain is set at the lowest acceptable 
value. Using the preferred Standard method of analysis, the uncompensated 
system is analyzed. The minimum gain permissible is that which provides 
acceptable steady state error. In most uncompensated systems this gain 
exceeds the critical gain of the system and the system becomes unstable. 
Proper compensation is necessary to make the system stable maintaining the 
required gain. Many different methods of compensation are available depend- 
ing on situations, but whenever possible, a phase lag compensation is recom- 
mended because the compensator is also acting as a smoothing filter. The 
required band width is narrow in this system because the parameter drift 


frequency is generally very low and it is effective in suppressing noise. 
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5.5 Example 

Some adaptive control systems are designed analytically and simulat- 
ed on computer. The actual performance of the system is compared with the 
one predicted by analytical method, Because of difficulty in simulation 
on an analog computer, the first example uses mild compensation and was 
simulated on analog computer. In the last example an extreme case of 
design is made and it is simulated on digital computer. 

(a) Example 1 
Design an adaptive control loop to maintain optimum for the 


following plant in which the pole is drifting with environment change. 


81 
G (s) = 3+9) 


Input to the plant; R = 4sin1,884t 
Specification; Steady state error = 10% of average velocity of the para- 
meter drift. 
Maximum overshoot = less than 20% 
Hunting loss = 10% of optimum parameter value 
Expected maximum frequency of parameter drift 
= 0.1 cps. 
Design steps 
(1) Gain requirement 


Steady state error = Wa 


K 
у 


where; A), = average velocity of parameter drift 
К, = static error coefficient 
In this case K, = 10 


(2) Selection of perturbation 
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Parameter perturbation signal; ES = С, sin Y t 
1 
Assume 1 v of parameter adjusting signal corresponds to 10% of optimum para- 


meter value. Therefore, 2 =lv 


The expected maximum frequency of parameter drift is given as 0 n = 0.1срѕ. 
Сһооѕе w, = 5 cps. which gives Ww: e = 1 : 50 and gives also no sharp 
overshoot in correction signal as seen in simulation study. Thus, E = 1 sin 
31.46 
(3) Establishing a linearized block diagram of the system 
Referring to the Fig. 4 the complete block diagram of the system, 
find each block. 


1) G(s) is generally, from 3.1 


m 
Gos) 7 1+0) 5 
m 
The measured phase characteristic of G (s) is the same as the phase character- 
istic of the plant (Refer to the section 4.1(c)(2) ). Therefore, the denom- 
inator of the G(s) is the same as the plant transfer function denominator, 


thus K 


m 
G8) = $13 0.1115) 


К 
m 


|6 (5) Se EE EE 
اھ‎ (jw ,) (1+0.111jw_) 


II 
> 
II 


0.0785 rad/sec. 


A was measured as 0.086 at О 
K is calculated as 0.122, thus 


0.122 
C OS) 7 $40.111$) 


2) M, 


The linearized M, is obtained from equation (5-3) 


Gy = 0 
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The a, is assumed to be 1. The attenuation factor of the perturbation 
signal B was measured. The value of 8 is different depending on the 
point of the error output signal wave at which it is measured as shown in 


Fig. 10(a). 


0.0087 


Maximum à: 


Minimum A - 0.00199 
Average i = 0.0054 


Take the maximum E for the first trial design, С, = 1 v, thus 


1 


Смт = 0.0174 


The value of ж: affects the overall loop gain which will be adjusted 


properly at the final design, 


"E 
3) с, (8) => 1 
S + (0,5 + 
1 1 
Linearized r (s) is, from equation (5-5), 
A 
— E er S 


A = A, = attenuation factor of G) (s) at W., which is measured as 6,6 


3 1 


А8 = : : = о 
Ze. (j ¿O + j 00) = 2.67 


tan 2.67" Тед $ 


S= j@_ = j0.628 
ран 2.67” 


"След. ~ 0.628 ہت‎ 
Thus, 
6.6 
Gea. 08) = 11000005: 
4) M, 
1 BC, 


From equation (5-8), Guo = 3 
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B is the attenuation factor of reference signal phase lag filter G, (s). 
In this example G, (s) is not used because the G (s) gives the carrier phase 
shift of 160° which is closer to 180° phase shift. Therefore, B = 1 


C, = 1 v. Thus, = 0.5 


1 бм2 
5) G, (s) will be a compensator in this sytem and is not decided yet. 
6) The integrator G, (s) is chosen to have following form because of diffi- 


culty in the simulation circuit as explained in section 4.2. 
K 
4 


9448) = 1+ 105 
(4) The linearized block diagram of the system becomes as Fig. 35(a) and 
which is reduced to Figure (b). Excluding G,(s) from the system and put- 
t ing ce = 10, the stability of the system is investigated. To obtain K, = 
10, K, should be 1430. The open loop transfer function of uncompensated 


system becomes, 


СН PU = 
о 5(1 + 105) (1 + 0.1113) (а + 0.06035) 


The Bode diagram of the system is plotted as Fig. 36. The phase margin is 


(5-9) 


zero and the system is unstable. 
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G (s) M, MO 









0.122 6.6 





$ (1+0.111$) 1+0.06035 







(а) 






0.007K, 


$ (1+105) (1+0.1115)(1+.0605$) 


(b) 


Fig. 35 The linearized block diagram of example 1. 
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The G (s) also can be written in standard transient response form as, 


Aia a E 
G, (9) 7 stg 0.1 (s « 906 + 16.6) a 
The root locus is also plotted as Fig. 37. 
The roots are; 
r, = 0.0336 - 10.9874 
М2 = -9.2322 
г, = -16.5349 (5-11) 


Two roots are located in the right half of the S plane. Thus the root 
locus plot also shows the instability of the system. 
The characteristic equation of the closed loop system is, 


g^ * 25.753 + 151.9874 14.98 + K=0 (5-11) 


By Routh's criterion, the critical gain of the system is found as К. = 87.9 
The required gain of the system is 149 and exceeds the critical gain. 

A compensation is necessary to make the system stable. 

(5) Compensation 


A lag compensation is adopted, the compensator is 


l + aTS 


G, (s) = TTT )5-12( 


Compensation 1; 

To get 5° phase margin in the final compensated system, a cross-over fre- 
quency а) = 0.5 (phase margin is 7.5° on Bode diagram) is chosen in Fig. 
36. 


aT is calculated from the following equation, 


1 + We (5-13) 
aT 10 
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1 = 0.5 ат = 20 (5-14) 


le (30) | = 12 db on Bode diagram Fig. 36 


This must be attenuated, therefore, 

20 log a = -12 а = 
Substitute this value in equation (5-14) and get T = 79.6. Thus, 
(5-15) 
This is the required compensator to have 5° phase margin. 
(6) Investigation of stability of the compensated system. 


Putting this G, (s) into the original system and you get the open loop 


transfer function of the compensated system as, 


O 101+05 (5-16) 
oc ~ §(1+79.6S) (1+10S) (1+0. 1115S) (1+0. 0603S) 
The Bode diagram of the compensated system is plotted in Fig. 36. Now, 
the system has 5° phase margin and expected to be stable. 
The open loop transfer function of the compensated system is rewritten 


in standard transient response form as, 


И 37.45 (5+0.05 (5-17) 
oc `°? = $5(5+0.0126) ($+0.1) (5+9) ($+16.6) 


The characteristic equation of the closed loop system is 


57425. 7138 14152. 224S°416.813S8-+37.638S+1.873 = 0 (5-18) 


The roots are 


-0.009526 + 10.4974 


м 
1 


-0.009526 - j0.4974 


rt 
tl 
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All roots are in the left half of the S plane. 
Appioximated +f of the compensated system is found as 0.01916 


M for this J = 1.9417 


(7) Simulation of the compensated system on analog computer 

The compensated system is simulated on analog computer. The actual 
response of the system is recorded as Fig. 38 where (a) is the parameter 
drift being corrected and (b) is the corresponding correction signal. The 
response of the system is different depending on the direction of drift of 
parameter as explained in Section 4,1(b). We take the worst case of response, 
that is, the response of the system when the parameter drifts to the direc- 
tion of decreasing system گج‎ . From the recorded result of the system 
response, following data are obtained, 


Steady state error = 16,51 


c - approximately 2 


(8) Comparison of theoretical and actual data of response 





theoretical actual 
steady state 
error 10% 16.57 
1.9417 2.0 
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(b) Corresponding correction signal 


(jt Parameter drift being corrected 


Fig. 38 The response of the system, Example 1. 
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(b) Example 1.1 

A different compensator is designed for the example 1 and the compensat- 
ed system is simulated on analog computer. The response of the actual system 
is also compared with the response predicted analytically. 
Compensation 2 ; 
(1) To get 10° phase margin in the final compensated system, Ww, = 0.3 
(phase margin is 17° on Bode diagram) is chosen from uncompensated system 
Bode diagram. 

Using equation (5-13), aT = 33.3 

G (jw) | = 20.5 db (From Bode diagram of uncompensated system) 


20 log a = -20.5 а = 106 


Т = 353 


Thus, Tu SENE (5-19) 
Go(s) = 353 5 


This is the required compensator to have 10? phase margin in the final 
compensated system. 

(2) Investigation of stability of the compensated system 

Putting the above G, (s) into the original system and get the open loop 
transfer function of the compensated system as, 


10(1+33. 35 š 
бос (8) = 5 +3535) (+103) (1+0.1115) @+0.06035) (5-20) 


The Bode diagram of the compensated system is plotted in Fig. 39. Now, 
the system has 10° phase margin. 
The open loop transfer function is rewritten in standard transient re- 


sponse form as, 


14.1 (5-0.03 
= Р 5-21 
Goc lS) = S(S+0.00284)(S+0.1)(S-9)(S+16.6) ( ) 
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The characteristic equation of the closed system is 


5 4 
5 +25, 702845 +151.97353+15.32852+14,14235+0.423 = 0 (5-22) 


The roots are 


гү = -0.02784 + 10.3025 
r, = -0.02784 - 10.3025 
r. = -0.03062 
r, = -9,0157 
r = -16.6008 


All roots are in the left half of the S plane. 
Approximate Y of the compensated system is 0.0915 
Me for this f = 1.758 
(3) Simulation of the compensated system on analog computer 

The compensated system is simulated on analog computer. The result 
is shown in Fig. 40, where (a) is the parameter drift being corrected, 
(b) is the corresponding correction signal, (c) is the demodulated output 
and (d) is the error between the plant output and model output. 
From the results: 

Steady state error = 11.17 


pt 
(4) Comparison of theoretical and actual data of response 





theoretical actual 
steady state 
error 10% ТАТТУ 
М 1.758 2.0 


pt 





(5) Simulation of compensated system on digital computer 
Because of difficulty in the simulation of the required compensator 


on analog computer, an approximate simulation of the required compensator 
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Fig. 40; Example l.l response of compensated system 
(a) parameter drift being corrected 
(b).Correction signal 

k Demodulated he 

d) Error E 
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was used in the simulation on analog computer in (3). 

The same system is simulated on a digital computer and the result 
is compared. Refer to the Appendix I for the FORTRAN program of the 
simulation. Fig. 41 is the graph of correction signal plotted from the 
numerical data output of the computer. If the Fig. 41 is compared with 
Fig. 40(b) there is a difference in signal shape. This is due to the ap- 
proximate simulation of the compensator in Fig. 40. 
(c) Example 1.2 

To get more phase margin in the system of Example 1, another compensa- 
tor is designed and the compensated system is simulated on a digital compu- 
per. 
Compensation 3: 
(1) To get 40° phase margin in the final compensated system, @ = 0.1 
(phase margin is 45° on Bode diagram) is chosen from uncompensated system 
Bode diagram. Calculated in the same way as in example 1, the required com- 
pensator is found as 


1 + 1005 


625) = 1 108 


(2) Investigation of stability of the compensated system 
Putting this C, (s) into the original system and get the open loop 


transfer function of the compensated system as, 


10(1+1005) 
Goes? = 4 
S(1+10 S)(1+10S)(1+0.111S)(1+0,0603) 








The Bode diagram of the compensated system is plotted in Fig. 42. Now, the 


system has 40° phase margin. 


The open loop transfer function is rewritten in standard transient response 


form as 
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КИ (5): - 1.49(S+0.01) 
oc Ë -4 
S(S+10 )(S+0.1)(S+9) (S+16.6) 


The characteristic equation of the closed system is 


3 


5 4 
57%25.70015 %151.902575 +14.9151952+1.491495+0.0149 = 0 


The roots are 


ri = -0.04353 + 10.08431 
r, = -0.04353 -j0.08431 
г, = -0.01108 

т, = -8.9951 

ro = -16.6069 


All roots are in the left half of the S plane. 
From the pair of complex roots, ¡$ of the system is found as 0.459 and 
the m for this value of Jf is 1.195. But this is no more true because 
these complex roots are no more dominant in the system. The dominant root 
is г. and real. Accordingly the expected response of the system should be 
an overdamped one and the time constant is large. 
(3) The simulation of the compensated system on a digital computer 
The compensated system is simulated on a digital computer. The FORTRAN 
programming for the simulation on a digital computer is shown in Appendix I. 
In calculating K, for the simulation, considering an overdamped re- 
зропзе, Аз = 0.00199 which is the minimum value of да is used and the K, 
1$ found as 6240 to get К, = 10. The result of the simulation is shown in 
Fig. 43 which shows the pole drift being corrected and the corresponding 
correction signal. The Fig. 44 shows the pole drift being corrected and 
the error. From the results, the steady state error can not be found ex- 
actly because the signal does not reach the steady state in 40 seconds used 
in the computer simulation. 
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(4) Comparison of theoretical and actual rssponse 
It is predicted theoretically that the coxnpensatec system response is 
overdamped and the time constant is large, The actual response shows the 


same pattern of response as predicted, The thecretical steady state error 


к? 


15 10%. The actual steady state error could not 
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to reach close to the theoretical one. 
(3) Example 2 

To test the validity of the analytical methed im an extreme case of 
design, the steady state error in erample 1 is changed es 1% end rigercas 
compensation is used, All other specifications are the came 2s example 1. 
(1) The required K, = 100 
(2) The open loop transfer function of oncompenseteó system is 


6 (з) ЖА к к т а EA 002 92-6 (5-23) 
© = <(1 + 105) Ц + 0.1115) ( + 0.06035) 





The Bode diagram of uncompensated system is plotted as Fig. 45, The phase 
margin is -2¢° and the system is unstable, 


Тһе С o 8) is also rewritten in standaré transient resgoense form as, 














— — (5-24) 


9)5 + 15.6; 


ER 


G (5) * scs 77 6.1308 + 


ә 


The characteristic equation of closed loop system l 


e 


6 9 : m A €^ 
S 4 25,73 + 151.98 + 14.95 + 


p^ 
Y ~ 
мә 
с> 
fi 

= > 


The roots are 


ту = 0.6049 + j 2,8461 
г, = 0,639 j 2.8461 
D 211.209 





By Routh's criterion, the critical gain is found as K. = 87.9. 
The required gain 1490 exceeds the critical gain. 
(3) Compensation 

A lag compensation is adopted. To get 40° phase margin in the final 
compensated system, cross-over frequency A = 0.1 (phase margin is 45? cn 
Bode diagram) is chosen from uncompensated system Bode diagram, Fig. 45. 
Calculated in the same way as in example one the required compensator is 
found as, 

1 + 100S 


G@,(g8) = —— — (5-26) 
$ 1+ 102$ 


(4) Investigation of stability of the compensated system 
Putting this G, (s) into the original system and get the open loop 


transfer function of the compensated system as, 


aaO OOS [_ — 
100(1+1005 (5-27) 


Ж) E 5 
5(1+10 5) (1+105) (1+0.1115) (1+0.06035) 


The Bode diagram of the compensated system is plotted in Fig. 45. Now, 
the phase margin is 39". The open loop transfer function is rewritten ir 


standard transient response form as, 


(5-28) 





G (C) = 515507006 


` 


The characteristic equatiom of closed system is 


59,25,7000157%151.90025752%14,9015157%1,4901295%0.0149 zt (5-29) 


The roots are 


г. = -0.04348 + 3 0.08428 
г) = -0.04348 - 10.08428 
г. = -0.01109 

г, = -8.995 


E re = -16.6069 A 
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All roots are in the left half of the S plane and are the same values 
as the roots in the Example 1.2. Thus the same response as in the example 
1.2 is expected, that is, the response of the compensated system is over- 
damped and the time constant is large. Furthermore, a smaller steady state 
error is expected theoretically. 

(5) Simulation of the compensated system on a digital computer 

The compensated system is simulated on a digital computer. The FORTRAN 
programming for the simulation on a digital computer is shown in Appendix 1. 
In calculating K, for the simulation, ES = 0.00199 15$ used and the K, is 
found as 62400 to get K = 100. Considering the overdamped response, K, 
is increased by 35% of the calculated value to 84000 in order to see the 
effect of an increased gain. The system still maintains stability and the 
rise time is reduced a little compared to the response of example 1.2. The 
result of the simulation is shown in Fig. 46 which shows the pole drift 
being corrected and the corresponding correction signal. From the result, 
che steady state error can not be found because the signal does not reach 
the steady state yet. 

(5) Comparison of theoretical and actual response 

it is predicted theoretically that the compensated system response is 
overdamped and the time constant is very large. The actual response shows 
the same pattern as predicted and also shows that the 35% more loop gain 
than the designed loop gain does not hamper the stability of the system at 
211, which tells that the system is in very stable condition as designed. 
The theoretical steady state error is 1%. The actual steady state error 
could not be measured but the deviation of the steady state error is ex- 


pected to be not so large. 
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6. Conclusions. 
(a) Analysis of the system. 

The open loop signal analysis showed clearly the feasibility of 
the adaptive scheme in this system, 

The signal analysis leads to the obtaining of a correct equivalent 
signal transfer block for the component in the system. 

The simulation study of the successive adaptive circuit, from the 
simplest to the complete circuit, showed the function of each component and 
the effect of component change in the circuit. 

It is confirmed that only the parameter drift signal component is 
acting through the closed adaptive loop; accordingly only the parameter 
drift signal need to be considered in the closed loop analysis. 

The simulation study showed that the analytical explanation of the 
system operations such as modulation, demodulation and filtering, are in 
close agreement with the actual operation of the system. 

(b) Design of the system, 

The obtaining of a correct error measurement transfer function 
G(s) and the attenuation factor of perturbation signal through the plant 
ا‎ is most important in the analytical design because these affect the 
accuracy of the analytical design the greatest. The value of E differs 
depending on the point of the error output wave at which it is measured. 
When the compensated system is expected to be underdamped, high value of 
A is used in the simulation and when the compensated system is expected to 
be overdamped, low value of En is used in the simulation for the first trial 
design. 

All the components are properly linearized and an analytical design is 


attempted. The predicted performance of the analytical design is in close 
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agreement with the actual system performance as verified in examples. 

In the case of a compensated system which has a pair of dominant 
complex roots (Example 1 and 1.1) the steady state error deviation is 
within 10% and Us deviation is within 15% between the predicted and actual 
value. 

In the case of a compensated system in which the pair of complex roots 
are not dominant and a real root is dominant, the comparison could not be 
made numerically, but the actual time performance had the same pattern as 
predicted. 

This analytical method of design simplifies the design considerably 
and is useful in obtaining a preliminary design which will be refined by 
an experimental verification. 

It is further suggested that this linearized technique is applied to 
an analysis and design of multi-loop sinusoidal perturbation adaptive con- 
trol system with some modifications for the intercoupling of adaptive loops, 


which should be further investigated. 
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APPENDIX I 
FORTRAN PROGRAMMING FOR THE SIMULATION OF 
A SINUSOIDAL PARAMETER PERTURBATION 
ADAPTIVE CONTROL SYSTEM 

To program the simulation of the system, the original block diagram 
of the system is rearranged to a proper form for programming as shown in 
Fig. 47 which is the rearranged block diagram of the complete adaptive 
control circuit, Fig. 32. Then, a simultaneous differential equations are 
set up in the form shown at the end of the Program Adapt l (page 109 ). 
The Program Adapt 1 in the following two pages, is the FORTRAN program for 
the simulation of the compensated system in Example 1.1. The input data 


and the graph control data are shown in proper form in page 110. A 


sample numerical output is shown on page 111. 
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